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ABSTRACT 



Aims. We report on X-ray observations of the supernova remnant 0509-67.5 in the Large Magellanic Cloud with XMM- 
Newton X-ray observatory. We use the imaging spectroscopy (EPIC) and Reflective Grating Spectrometer (RGS) data 
to investigate properties of the remnant and its environment. 

Methods. The X-ray spectra were analyzed with SPEX software package. In addition to this we performed a numerical 
hydrodynamic simulation of the remnant. 

Results. The EPIC data show prominent Fe K line emission, but the deduced overall amount of iron in the shocked 
ejecta is low. The data also show that the remnant has an asymmetric ejecta structure: the bright southwest region of the 
remnant shows an overabundance of metals. The analysis of the RGS spectrum shows that the remnant has a high lines 
velocity broadening of 5000 km/s. We found a hydrodynamical model for the remnant with basic hydrodynamical and 
spectral parameters similar to the observed ones. 

Conclusions. The data analysis show that the reverse shock just recently reached iron layers of the ejecta. The bright- 
ness enhancement in the southwest region could be a sign of an asymmetric explosion or it could be the result of a 
density enhancement of the interstellar medium. We constructed numerical models which are in good agreement with 
the observations, with circumstellar density of 3 x 10"^^^ g/cm^, age of ~ 400 years, velocities of ~ 5000 km/s and an 
electron to ion temperature ratio of 10"^. 



Key words. X-rays: individuals: SNR0509-67.5 - 
nants — method: data analysis — hydrodynamics 



1. Introduction 



The X-ray emission from young supernova remnants 
(SNRs) provide a vital source of information about the 
properties of shocked, rarefied plasma, properties of the 
ambient medium around the remnants, and also about 
supernova explosion models. With the current genera- 
tion of X-ray telescopes, XMM-Newton, Chandra, and 
Suzaku, one can combine now the imaging and spectral 
information to study the extended emission from SNRs. 
Moreover, in the case of Chandra and XMM-Newton, dis- 
persive spectrometers allows one to obtain high resolution 
spectra, providing more constraints on the plasma param- 
eters, and on the dynamics of the plasma, through Doppler 
shifts and broadening. The XMM-Newton Reflection 
Grating Spectrometers (RGS) have the advantage that the 
spectral quality is less detonated by the spatial extent of 
the object, at least for sources with an angular extent of 
< 1 '. 



ISM: individuals objects: SNR0509-67.5 — ISM: supernova rem- 



For that reason the Large Magellanic Cloud (LMC) 
and the Small Magellanic Cloud (SMC) SNRs are par- 
ticular well suited, since the relative proximity of the 
LMC means that many SNRs are bright enough to ob- 
tain high signal to noise spectra with both the CCD and 
grating spectrometers. On the other hand, the distance of 
the LMC, 50 kpc, is such that young SNRs have an ex- 
tent less than 1 ', making them excellent targets for the 
XMM-Newton RGS. The XMM-Newton CCDs, with a 
resolution of ~ 6 ", still allows one to study spatial varia- 
tions, although Chandra obviously provides more spatial 
details. In addition, studying SMC/LMC remnants has the 
advantage that the interstellar absorption is on average 
lower than for most of the Galactic remnants, and their 
distance is known. 

XMM-Newton observations of several LMC/SMC 
remn ants have been pubUsh ed: N132D j Behar et aTl 
20011). IE 010 2.2-7219 ("Rasmussen e t al.l 12001 



Sasaki et al.l |2006|) . B0540-69.3 (ivan der Hevden et al.l 
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1200 lb . N103B (^van der Hevden et alj 



SMC remnants (van der Hevden et al. 



2002) 



2004 



I, and several 2. The XMM-Newton data, reduction and 
analysis 



Here we report on the X-ray spectrum of the 
SNR 0509-67.5 in the Large Magellanic Cloud. This 
object was alr eady examined by several auth ors, e.g.. 



IWarren & Hughes (2004); Vink (2006); Badenes et al. 




iugi 



Rest et al.1 (I2008h; iGhavamian et al. 



Rest et al 



Badenes et al.l (I2007h ; 



iHendrick & Revnoldsl ( 1200 Ih . The remnant has a spher- 



ical shape and is somewhat similar in age and size to 
Tycho's SNR, which is also a Type la SNR. SNR 0509- 
67.5 has an angular dia meter of 25", correspon ding 
to the radius of 3.6 pc dWarren & HughesI |2004 . A 
recent analysi s by Rest et al. (l2005b . based on light echo 
detection, and iBadenes et al.l ( |2008|) . based on modeling 
of the X-ray observation, sugge st the age of the remnan t 
to be ~ 400 years. Estimates bv IGhavamian et al] ( |2007|) 
gives the values of 295 - 585 years, consistent with the 
rough estimate ~ 500 years by IVinkI (l2006h . based on a 
preliminary analysis of the XMM-Newton data. 

A detailed analysis of the X-ray emiss i on of 
SNR 0509-67.5 was made by I Warren & HughesI (|2004 
and based on Chandra X-ray data. Their study revealed 
that the remnant is rich in silicon, sulfur and iron, the bulk 
of the continuum emission has a non-thermal origin, and 
an estimate the circumstellar medium (CSM) density is 
of ~ 0.05 cm"^. The preferred model for the explosion 
which produced the remnant — is delayed detonation. 



Recently, iBadenes et al.l (l2008l) reported their study 
of the remnant's X-ray spectrum observed by Chandra 
and XMM-Newton observations. They performed hydro- 
dynamical and X-ray spectral calculations of the remnant 
and concluded that it was an energetic, delayed detona- 
tion explosion of 1.4 x 10^' ergs with a nickel mass of 
O.97M0. In their numerical simulations the circumstellar 
density was set to ncsM - 0.4 cm"^ and the ratio of elec- 
tron to ion temperatures to 0.02. 

For illustrative purposes we show in FiglT] some of 
the Chandra X-ray images of the remnant in three differ- 
ent bands. The overall structure of the SNR shows a large 
scale inhomogeneity in the shell. The Southwest (SW) 
part of the remnant is somewhat brighter compared to the 
average brightness. 

In the current study we performed an analysis of the 
remnant employing XM M-Newton X-ray ob servations 
and SPEX fitting software dKaastra et al.lll996l) . We also 
performed hydrodynamical simulations of the SNR 0509- 
67.5, employing the hydrocode supremna dSorokina et al.l 
|2004 to compare numerical predictions for the remnant 
with the observations. 

The paper is organized as follows. First we describe 
briefly the XMM-Newton data in ^ then we discuss the 
X-ray spectrum in ^ Numerical models for the SNR are 
presented in ^ We discuss our results and outline them 
in^ 



SNR 0509-67.5 was observed with the XMM-Newton X- 
ray Observatory (obs ID 01 1 1 130201) for 35.9 ks on July 
4, 2000. The observatory contains three X-ray telescopes. 
Behind all of them are CCD detectors, collectively called 
the European Photon Imaging Camera (EPIC). Two of 
them ar e MPS type CCD detectors called MOSl and 
MOS 2 dTurneret al.1 1200 lb . the other is of the pn-CCD 
type dStriider et al.ll200l1) . The telescopes equipped with 
MOS CCD detectors in the focal plane contain also the 
Reflective Grating Spectrometers (RGS), consisting of 
the Reflective Grating Arrays, dispersing about 50% o f 
the X-rays to two CCD arrays dden Herder et alj|200lh . 
XMM-Newton observations therefore provide simultane- 
ous data from all five X-ray instruments: the two MOS- 
CCDs, the pn-CCD, and the two RGS instruments. 

In this paper we concentrate on the data obtained with 
the MOS detectors and the RGS. Although the EPIC- 
MOS have a somewhat lower sensitivity than the EPIC-pn 
instrument, they have better spectral resolution, which is 
important for line rich sources, such as SNRs. 

Note that IBadenes et al.l d2008l) found an inconsis- 
tency between EPIC MOS 1,2 and PN spectra. They in- 
vestigated Si Ka line centroid location in the datasets of 
these devices and found out that the line centroids of 
MOS and PN detectors are shifted with respect to each 
other. The authors rejected MOSl and MOS2 data and 
concentrated only on the Chandra and EPIC-pn spectra. 
We confirm this inconsistency: analysis of the Si Ka cen- 
troid of MOS 1,2 and PN data reveals that it is located at 
1.8503 + 0.0012 keV and 1.8298 + 0.0015 keV respec- 
tively (the errors are Icr x^)- For the Fe Ka line we found 
6.495 + 0.048 keV and 6.442 + 0.025 keV. The discrep- 
ancies in the lines locations are about 1%. Nevertheless, 
there are indications, that PN instrument has small gain 
problems and sometimes a shifting of the PN energy 
grid is required (private communicatio ns with N.Werner; 
I Werner et aT]|2006l; lde"Plaa et al.ll2004 . Moreover, in our 
study, we were unable to find a fitting model with sensible 
parameters of the emitting plasma for the PN data, mean- 
while a spectral model for the MOS data gives reasonable 
and satisfactory plasma parameters. 

The EPIC data was slightly affected by background 
flaring ("soft protons"). As a result we cut out 3.4 ks of 
the total observation. The MOSl observation was done 
in small window mode, whereas the MOS2 data were 
made in full frame mode. For MOSl, therefore, we could 
not select a region to extract background spectra. For the 
MOS2 this was possible, but we found that the back- 
ground correction was too small, and the account for it 
did not change the fitted parameters. For that reason we 
did not include background subtraction for our final anal- 
ysis. 

The RGS is a slitless spectrometer For an extended 
source this means that the spectrum is smeared by the 
image of the source itself. For SMC/LMC remnants the 
smearing is modest, but present, and it gives rise to 
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Fig. 1. Chandra images of the SNR 0509-67.5 in different 
X-ray ranges. Fr om left to right: 0.45-1.7 5 keV, 1.8-2.5 
keV, 2.5-6.0 keV (IWarren & Hughesll2004 . 



a change in the line spread function. For our analysis 
we incorporated this effect into the response matrix by 
convoluting the standard (point source) response matrix 
with the emissivity profile of the SNR, as obtained from 
archival Chandra observations. This procedure was also 
for the RGS data of SN 1006 (Vinket al. 2003). As we 
shall show later, however, in the case of 0509-67.5 the 
Doppler line broadening is much larger than broadening 
due to the spatial extent of the remnant. 

Apart from adapting the RGS response matrix, all re- 
duction for both MOS and RGS data was done using 
the standard XMM-Newton software package SAS ver- 
sion 7.1.0. 



3. Spectral models 

The X-ray spectra of the SNR were fitted with a non- 
equilibrium ionization ( NEI) model, which is a part of 
the SPEX fitting package dKaastra et al.lll996l) . The pack- 
age does not include a plane shock model (which takes 
into account temperature and ionization timescale gra- 
dients), but it contains the most complete and up-to- 
date set of atomic data. The use of a single ionization 
timescale NEI model is justified in some cases, when 
one nee ds to get rough estimates o f basic SNR proper- 
ties (e.g. lCassam-Chenai et al.ll2004h . The emission mea- 
sure nnrieV in such a model can be used to estimate the 
circumstellar medium density. Since continuum emission 
stems mostly from the shocked CSM, the value of hh can 
give us directly an estimate of the unshocked environment 
density (depending on the assumed equation of state). 

Combined fitting of the spectra from RGS, EPIC: 
MOSl, MOS2 devices is presented in Fig|2] 



3.1. EPIC data 

3.1 .1 . Spectra parameters 

In order to adequately fit the spectra we needed at least 
two thermal components: one NEI component is for the 
bulk of the X-ray line emission, and another one is 
to fit the Fe-K emission around 6.5 keV (see section 
13.1.21 for details on this component). We also investi- 
gated ajn£del_wath_an additional power law continuum, 
since lWarren & HughesI ( |2004|) reported indirect evidence 
for non-thermal continuum, probably synchrotron emis- 
sion from > TeV electrons. The plasma parameters for the 



two models are indicated in Table[T]as (th) and (th+pow). 
In addition, we extracted also the spectrum from the en- 
hanced SW region of the remnant. In general, the pa- 
ra meters are in agreement with the the ones obtained 
by IWarren & HughesI (|2004|) . The enhanced SW region 
shows higher ionization timescale, probably reflecting a 
higher density. 

Using the value of the emission measure, we can es- 
timate the density of the CSM. Hydrogen emission mea- 
sure (Table [U MOS+RGS) is EM = «,«« V = 1.5 x 



10 cm and the remnant's size is /? - 3.6 pc ^ 
lO'^ cm. Assuming that the emitting shell has a width of 
A/? = /J/12, so that the volume Vx = V,or/4 = nR^ll, ^ 
10^^ cm-', we obtain «e«/f ^15 cm"^. Furthermore, as- 
suming that all the matter is ionized «h - 4 cm"-', we 
obtain that CSM number density for the remnant should 
be «csM - 1 cm"^ for the neutral circumstellar environ- 
ment, which appears to be the case for SNR 0509-67.5 
dGhavamian et al.ll2007h . Here we assumed shock com- 
pression ratio to be 4 for non-relativistic matter with adi- 
abatic index y - 5/3. However, cosmic ray (CR) acceler- 
ation may play an important role in SNR dynamics. For a 
cosmic ray dominated shock the adiabatic approach that 
of a relativistic gas (y - 4/3) gives rise to a compres- 
sion factor of 7. Cosmic ray e scape may even increase th e 
compression ratio beyond 7 dBerezhko & Ellisonlll999l) . 
Tycho's SNR provides s ome observ ational evidence for 
high compression factors (IWarren et al. 2005) . Therefore, 
in the case of 0509-67.5, the pre-shock density could be 
as low as «csm ^ 0.6 cm"-'. 

Also, two additional sources of uncertainty concern- 
ing the pre-shock density are the presence or absence 
of a non-thermal continuum, and chemical composition. 
In particular, our above estimate assumes that most of 
the continuum is bremsstrahlung from hydrogen-electron 
collisions, whereas it is quite likely that parts of the 
remnant consist of pure metal plasmas, in which case 
bremsstrahlung comes from metal ions and the electrons 
they provide, in general this lowers the density estimate. 
Therefore, our estimate of «csm - 1 cm""*, should be con- 
sidered as an upper limit. 

Abundances of chemical elements for the entire 
SNR (thermal and nonthermal continuum) and enhanced 
SW part for the simple NEI model are presented in 
Fig [3l Abundances for some theoretical explos ion mod- 
els llwamoto et al.1 119991: IWooslev et"an l2007h are also 
shown. The accounting for the nonthermal continuum 
does not create a considerable difference in the abun- 
dances set, but naturally it makes the remnant more metal- 
rich. One can see from the figure that the SW region is 
more metal abundant compared to overall remnant values 
(2-3 times higher amount of metals). 

The SPEX fitting shows no evidence for interstellar 
absorption, which means that the value of Nh is low 
for this remnant. The c olumn density value found by 
IWarren & HughesI ( |2004|) is Nh ^ {5-9) x \(fi° cm^. 
Which is, indeed, the lo west value among all known LMC 
supernovae remnants dHughes et alj |I998|) . This could 
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Fig. 2. MOSl, MOS2 and RGS data fitted with the NEI 
model, the Fe K line is fitted separately: a NEI compo- 
nent with kT - 3.5 keV and Hgt =10' s/cm^. Vertical 
lines show the location of three Gaussians, which substi- 
tute missing in spex model Fe L lines (the energies are 
1.102, 1.125 and 1.226 keV). 



suggest that the remnant is located on the near side of 
the LMC. 

The power law model for the continuum gives a 
value for the nonthermal flux of Ffj^^y = (2.3 + 0.3) x 
lO""* ph/s/cm^/keV with power index of F = 3.5 + 0.1, 
which is steeper than expected for this kind of rem- 
nants, e.g. Tycho SNR (c . f. Tycho's SNR or C as A 
ICassam-Chenai et al.ll2006l: IVink & La ming"20 d3|). This 
value is close to the one obtained by iWarren & HughesI 
(|2004 : 3.25 + 0.18. The extrapolat i on of the non-fliermal 
radio flux JHendrick & Revnoldsl 1200 ll) to the X-ray 
range yields Fucev = 1.4 x 10"^ ph/s/cm^/keV. Since 
the X-ray synchrotron spectrum is likely to be affected by 
synchrotron losses, this flux should be considered as an 
upper limit. 



3.1.2. Fe-K emission 

Despite of the low amount of iron in the shocked ejecta 
deduced with single ionization timescale NEI model, the 
spectrum has a pronounced Fe K feature. Its centroid is 
located at 6.50 + 0. 05 keV, which may cor respond to 
the ions from Fe VII (iKaastra & Mewe|[l993h with riet - 
lO"* s/cm^ up to Fe XVII with n,t < 10'° s/cml We fit- 
ted the Fe K line separately with the temperature kT = 3.5 
keV and ionization timescale «ef = lO' s/cm"*. The ion- 
ization time is uncertain, but it is impossible to fit the Fe K 
emission with a component with net > 2xl0's/cm^, since 
then the spectrum would have exhibited much more Fe L 
emission than observed. Given the low n,,! it is likely that 
the Fe K emission is associated with pure Fe, recently 
shocked by the reverse shock. 

The separate fitting of the Fe K feature allows us to 
make a crude estimate of the amount of swept up iron 
layer of the supernova ejecta. Assuming that the pressure 




15 20 

Element number 

Fig. 3. Abundances of the entire remnant and its SW 
region, obtained in fitting of the EPIC (MOSl and 
MOS2) spectrum. A bundances for the classical W7 
dlwamoto et al.l 1 19991) and for the de layed-detonation 
model c050403m (IWooslev et al.l 120071) are also shown. 
Theoretical abundances are taken from numerical simu- 
lations (described below) include swept up (~ 0.7 M©) 
ejecta and shocked LMC circumstellar medium. All abun- 
dances are solar normalized. Errors are 1 cr rms. The sili- 
con abundance (with respect to solar) for the entire SNR 
is 5.0 and for the SW part of it is 15.7. 



and temperature do not change drastically throughout the 
shell, and thus the electron density n^ remains approx- 
imately constant, we can employ the following relation 
npeVpe ~ {EMFjEMx){nHVx). 

The value of the emission measure of the pure iron 
NEI component, /sMpe = n^npe Vpe - 4x lO^'* cm"-', basic 
NEI component, EMx - nenuVx 



10^^ cm ^, and also 



the density ne - 4 cm""* and the volume Vx - 10^^ cm"' 
of the emitting shell, yield us the estimate of the swept 
up iron in the remnant of Mpe = 56 my npeVpe ~ 0.05 M© 
(with mu = 1.66 X lO"^^ g, Mq = 2 x 10^^ g). This is 
less than 10% of the iron amount in a typic al la SN model 
JNomoto et alJll984l;llwamoto et al.|[T999l) . 



3.1 .3. Tine line emission deficit around 1 .2 keV 

The EPIC (MOSl,2) models do not fit an excess of emis- 
sion around 1.2 keV This feature is also present in the 
EPIC-pn spectra, and is weakly visible in the RGS spec- 
tra. Therefore, this excess is likely to be real and not 
caused by calibration errors or background features. It 
is possible to fit this region of the spectrum, but only if 
we allow for unrealistic overabundances of either Na (ap- 
proximately 5 times of Si) or Ni (4 times of iron). 
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Table 1. Parameters for the SNR 0509-67.5 (last column contains data for the southwest region), EPIC (MOSl,2) and 
RGS data. Errors are 1 sigma (;^f--distribution). th — basic (single ionization timescale) NEI model, th-i-pow — NEI 
with power-law continuum. 



Parameter 


th, MOS-hRGS 


th, MOS 


th+pow, MOS+RGS 


th-l-pow, MOS 


SW, th, MOS 


n.HH V X 10^^ cm-3 


1 15+0.12 


0.98-l« 


iS+ois 


95+0.11 


07+'"'5 


kT, keV 


4 Q 1+0.23 


5-02!°;i 


4 55+0.22 
^■-'•'-0.20 


4 98+o» 

^•'^''-0.30 


4 59+0.54 


n,fxlO'°, s/cm' 


1 41+0.03 

^■^^-0 03 


1 62+""^ 
i.oz„„, 


I 53+0.04 


1 60+""^ 

^•°"-0 0fi 


1 co+0.08 

^■-'^-oos 


cry,km/ s 


6030 ± 170 


5700 ± 220 


5050 + 180 


5350 ± 240 


6100 ±500 


X'ld.o.f. 


2.61 


1.42 


2.29 


1.37 


1.66 



In our view the possible cause is the uncertainty in the 
atomic data base of SPEX concerning the Fe-L emission. 
The Fe-L emission of SNR 0509-67.5 is dominated by Fe 
XVII emission. So it is likely that the excess is due to in- 
completeness of the atomic database concerning Fe XVII, 
in particular due to line emission from the « = 4 to « = 2 
level. This is corroborated by the fact that the series limit 
of Fe XVII is at 1.26 keV, which also roughly marks the 
end of the emission excess. Some « = 4 — > 2 lines are 
present in the SPEX code, but their emission strengths 
may be underestimated. This may not be surprising, since 
the code has been tested primarily on collisional equilib- 
rium sources, whereas in NEI n = 4 — > 2 may be en- 
hanced due to non-equilibrium effects, such as inner shell 
ionization of Fe XVlQ 

3.2. RGS data 

The high resolution RGS spectra allowed us to evaluate 
the velocity broadening of the X-ray emission lines and 
make an independent estimate of the CSM density. Thus, 
from the spectrum in energy range of 0.3 - 2.5 keV we 
find 0",, = 4900 + 420 km/s (which is the preferred value 
since it is based on the RGS only). This gives us an upper 
limit on the age of the SNR of t < /^/vshock - 540 + 50 yr 
(assuming Vshock = 4/3cr,.). 

The first column in Table [T] shows the best fit param- 
eters for the joint MOS-i-RGS analysis. However, we also 
fitted the RGS spectra separately with a single NEI model. 
In this case the fitted parameters are quite different, in 
particular the plasma temperature: kTg = 0.75^QQg keV 
and ionization timescale n^t - (1.28+^33) x 10"^ s/cm"'. 
The probable reason is that the RGS is mostly sensitive 
to the energies in 0.5-1 keV, so it misses the contribu- 
tion from lines that mostly emit at higher temperatures 
and the continuum. In particular the OVII lines around 
0.55 keV/22 A has much more statistical weight in the 



' We discussed this issue with Dr. Kaastra, one of the authors 
of the SPEX code, who agrees that this may explain the lack of 
a good fit around 1.2 keV. 



RGS spectra. The OVII line emission probably has con- 
tributions from a lower temperature plasma, that is not 
picked up by the MOS spectra. The discrepancy there- 
fore points to the presence of temperatures gradients in 
the SNR. 

Ideally, one would like to incorporate those gradients 
into a more complete model of the X-ray emission from 
SNR 0509-67.5. We made an attempt to obtain such a 
model by fitting a joint MOS-RGS model with 3 NEI 
components, one for the low temperature plasma, one for 
the bulk of the MOS spectrum, and one for the Fe-K line 
emission. However, it turned out to be difficult to come up 
with a unique, satisfying solutions, given the complexity 
of the parameter space. Moreover, it would still involve 
some arbitrary choices. For example, it is not clear what 
to do with the abundances; assume one set of abundances 
for the first two NEI components, but then we do not al- 
low for possible abundance gradients that are also likely 
to be present. So in short, we have reached here the limit 
of what is possible with multiple components. The RGS 
data, however, do show that the models in Table [T] are in- 
complete. 

One of the advantages of the high spectral resolution 
of the RGS is that it is possible to obtain an estimate of 
the nitrogen abundance. Nitrogen is not an ejecta product 
of Type la supernovae, so any nitrogen must come from 
the shocked CSM. Therefore, the nitrogen abundance, to- 
gether with the emission measure makes it possible to ob- 
tain an alternative estimate of the pre-shock CSM density, 
avoiding some of the confusion due to mixture of pure 
metal, ejecta components. 

The hydrogen emission measure of the RGS fit is 
of rieriH V - 4 X 10^^ cm"-' and emitting volume of 
Vx - 10" cm^ yield [n^nn]"^^^ - 400 cm^^ Taking into 
account best-fit RGS amount of nitrogen X^*^^ = 0.01 
and typical LMC abundance X'^^^ =0.4 (both in solar 
units), we derive for the circumstellar matter [«e«/f l*^^*^ = 
[n.riH]^'^^ (X^'^^/X^^^) = 10 cm-^ Thus for the shocked 
CSM (LMC abundances) «h = 3 cm"-', which yields 
ncsM = 0.4 - 0.8 cm"^ 
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Fig. 4. Top: the RGS spectra with best-fit NEI 
Bottom: the RGS spectra and the MOS-hRGS best- 
model. 
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4. Numerical models 

To simulate the evolution of the remnant and to com- 
pare it with the observed data, we employed the hy- 
drodyn amical code supremna , which is explained in de- 
tails in ISorokina et al.l (l2004h . The code assumes spher- 
ical symmetry, but it incorporates many relevant phys- 
ical processes, without which it is not possible to ac- 
curately predict peculiarities of SNR radiation, such as 
time-dependent ionization, inner-shell collisional ioniza- 
tion, possible difference in temperatures of electrons and 
ions, the influence of radiative losses, the account of 
electron thermal conduction and nonthermal particles. 
The calcu lations of the lines e mission was based on 
iGaetz & Sa lpeter (1983) and M ewe et al.l (11985). The in- 
tr oduction of inner -shell ionization processes is described 
in lKosenkol (120061) . 

To reproduce the observed remnant in numerical sim- 
ulations we considered various physical conditions of the 
CSM and also different explosion models: deflagration — 



W7 (iNomoto et al.l 19841) and some of the 'mildly-mixed" 
delayed-detonation models (Wooslev et al] |2007l) . 

Following the estimates of the CSM density derived 
from the spex spectral fitting, we surround the remnant 
with uniform medium of pcsM = 3 x 10"^^ gjcnr'. In 
this set up we found that at the age of f ~ 400 years the 
modeled remnant's radius reaches R ~ 3.6 pc, with typi- 
cal velocities of the plasma of v ~ 5000 km/s. These pa- 
rameters are in agreement with the observed ones of the 
SNR 0509-67.5. Modeled temperatures and X-ray spectra 
will be discussed in details below. 



4. 1 . Basic modeling 

Fir st, we used a "basic" approach, described 

in ISorokina et al.l ( |2004|) where electrons and ions 
are being heated partially due to artificial viscosity 
and degree of the electrons heating is controlled by a 
parameter qi, such that P,- = Pj{thermal) + qiQ and 
Pe - P eithermal) + (1 - qi)Q, where P^, Pi — electrons 
and ions pressure, Q — artificial viscosity (is the case 
when only collisional energy exchange takes place 
q, = I - nie/mi, m^, m, - electron and ion masses). Note 
that since q, controls the efficiency of electrons heating 
(their temperature may vary theoretically from 10 "^T, to 
Ti), this parameter has a strong inffuence on the X-ray 
spectrum behavior (line ffuxes with respect to free-free 
continuum). 

In this framework we have found that the best fit (in 
terms of x^/d.o.f.) for the observed spectrum gives the 
value of qi = 0.99, which yields the typical ratio of 
Te/Ti ^ 3 X 10"^ for the case of W7 explosion model. 
We also co nsidered a library of t hermonuclear explosion 
models by IWooslev et al.l (|2007|) . We employed several 
of them with various amounts of iron and intermediate 
mass elements (IMF). We have found that the explo- 
sion of a "mildly-mixed" delayed detonation model with 
£ = 1.4x 10^1 ergs. Mm = O. SMp, Mpe = QA M„, Mime = 
0.3Mo (c050403m model in IWooslev elall 12007, nota- 
tion) gives the best-fit to the observed FPIC MOS spec- 
trum of the SNR with qi = 0.9 and thus TJTi ^3x 10"^. 
These values seem to be in reasonable agreement with 
the relati on Te/Tj oc v~^ (vy — shock velocity) pro- 
posed by iGhavamian et aLJ ( l2000l) . Indeed, assuming for 
the SNR 0509-67.5 shoc k front speed of v ^ - 4 /3 v - 
6000 km/s, according to IGhavamian et al.l (l2000i Fig. 2) 
we should expect T^/Ti - 6x 10"^. 

The simulated (for these two explosion models) and 
the observed X-ray spectra are presented in Fig|6| Note 
that all of our models exhibit noticeable lack of ffux in 
the region of Fe L complex (compared to the observed 
spectrum). This is caused by some incompleteness of the 
atomic data set for Fe L emission employed in our pack- 
age, as verified by directly comparing a single NEI model 
as obtained by the supremna code and in spex. This issue 
will be addressed in future updates of the supremna code. 

At higher densities of the CSM, there is no the Fe L 
emission deficit, but in these models the expansion rate of 
the shell and the radius of the contact discontinuity are in 



Kosenko et al.: The X-ray Spectra Of The SNR 0509-67.5 



350 years 



400 years 




f^CSM 



CSM' 



Fig. 5. Maximum velocity of the shocked ejecta (v™?'^, 
upper row) and position of the contact discontinu- 
ity {RcD, lower row) depending on the value of 
CSM density for 350 (left column) and 400 years 
(right col umn). Explosion mo del — delayed detonation 
C050403 dWooslev et al.ll2007l) with £ = 1.4 x 10^' ergs. 
The observed values of the velocity and radius of the rem- 
nant are also outlined. 



disagreement with the observed values of the lines broad- 
ening and the remnant's size. So the value of CSM density 
was chosen to satisfy the observed velocity of the emitting 
ejecta, the age and the geometrical size (assuming that the 
observed edge of the remnant close corresponds to the lo- 
cation of the contact discontinuity) of the remnant for the 
considered explosion mechanisms. Figure |5] shows max- 
imum velocity of the shocked ejecta and position of the 
contact discontinuity depending on the assumed value of 
CSM density for the delayed detonation explosion model 
C050403 at the age of 350 (left panels) and 400 (right pan- 
els) years. The figure indicates that for an assumed age of 
400 years, the CSM density of 3 x 10"^^ g/cm^ gives 
a reasonable compromise between velocity and radius. 
An assumption that the remnant may be as young as 350 
years, allows the CSM to be denser than 5 x 10"^^ g/cm^, 
but in this case the synthetic X-ray spectra disagree with 
the observed ones (there are major discrepancies in Fe L 
emission and Si K and S K lines fluxes). More energetic 
(E > 1 .4 X 10^') progenitor could shift the allowable val- 
ues of the density to higher ranges. We checked several 
of them, which contain more ^6A^/, and which could be 
consistent with the optical sp ectrum of the sup ernova, as 
extracted from the light echo (iRest et al.ll2008h . However, 
these models contain too little Si and S to fit their emis- 
sion lines. 

X-ray spectra of the simulated and the observed 
(merged MOSl and MOS2, crosses) remnants. The the- 
oretical spectra are based on W7 (Te/Tj ^ 0.003, 
X^/d.o.f. ^ 27, dashed line) and delayed-detonation 
c050403m (rjr, ^ Q.Q3, x^/d.o.f. ^ 19, solid line) ex- 
plosion models, pcsM = 3 x 10"^^ g/cm^, age t - 400 
years. 




Fig. 6. X-ray spectra of the simulated and the observed 
(merged MOSl and MOS2, crosses) remnants. The the- 
oretical spectra are based on W7 {Te/Ti ^ 0.003, 
X^/d.o.f. ^ 27, dashed line) and delayed-detonation 
c050403m (TJTi ^ 0.03, x^/d.o.f. ^ 19, solid line) ex- 
plosion models. pcsM = 3 x 10"^^ g/cm"*, age t = 400 
years. 



The Table |2] shows a list of some basic properties of 
these models. Some of the numerical values in the table 
are upper limits or just typical for the corresponding case. 
Errors are presented where appropriate. 

4.2. Accounting for the cosmic rays 

We introduced an additional component representing the 
pressure and energy from (relativistic) cosmic ray parti- 
cles, in a way that could be easily implemented within the 
current framework of this code. This was done by altering 
(softening) the equation of state as follows; 



P = 
E = 



Pg + iar^ 

Eg + aT'^ 



(1) 



where Pg is the gas pressure. Eg is the gas internal en- 
ergy, T is the "temperature" and a is a free parameter 
mimicking the contribution of the relativistic particles. 
Effectively it means, that some part of the energy is trans- 
ferred into relativistic particles, which behave similar to a 
photon gas. 

Density profiles for the cases with and without 
this relativi stic correction for the delayed-detonation 
c050403m dWooslev et al.1 12007|) model are presented 
in Fig|7] Note that in the case of the relativistic alteration 
of the equation of state (gray solid line), the density jump 
at the forward shock is enhanced, compared to the basic 
set up. The density profile of the ejecta behind the reverse 
shock became affected as well. 
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Fig. 7. Density profiles (solid lines) against radius for the 
hydrodynamical models of age 400 years and for the de- 
layed detonation explosion mechanism c050403m. Upper 
panel; basic set up, lower panel: the set up with the rela- 
tivistic correction in the equation of state (Eq.[TJ. Also the 
distributions of Si (black dash-dotted line) and Fe (gray 
dashed line) in the ejecta are shown. 



Fig. 8. X-ray spectra of the simulated and the observed 
(merged MOSl and MOSl) remnants. The theoretical 
spectra are based on W7 (x^/d.o.f. ^ 21, dashed line) 
and delayed-detonation c050403m (x^/d.o.f. ^ 13, solid 
line) explosion models. pcsM = 3 x 10"^^ g/cm"*, age 
f = 400 years, EcR/Eg ^ 0.4. 



However, in reality the cosmic ray contribution to the 
internal energy may be different for the forward and re- 
verse shock. Recentlv, iHelder & Vink ( 2008) established, 
that in the case of the SNR Cas A, the reverse shock can 
accelerate cosmic ray particles up to energies of a few 
TeV. 

In this set up, a configuration with the cosmic-ray en- 
ergy density Eqr - ciT'^ of the order of 0.4£'^, produce a 
model which also gives a good fit to the observed spec- 
trum. The modeled and observed X-ray spectra for this 
approach are presented in Fig IS] 

I Chevalier! (Il983h considered self-similar solutions for 
a mixture of a fluid with 7 = 5/3 and a relativistic gas 
(i.e., cosmic rays) with y - A-jli. The effective adia- 
batic index was determined by a constant fraction of cos- 
mic rays in total pressure. Later lBlondin & EUisonl ( 1200 ll) 
used a constant value of effective j which was taken 
sometimes even lower than 4/3 to allow for a stronger 
compression ratio mimicking radiative shock waves. We 
note that our parametrization, being also simple in com- 
parison with a true cosmic r a ys model, is much richer 
than used bv lChevalieij ( ll983l) : lBlondin & EUisonl ( 1200 ll) 
in their purely hydrodynamic studies. We do not assume 
a constant fraction of cosmic rays or a constant value of 
effective y, thermal conduction and radiative losses are 
fully included into a hydrodynamical scheme, in contrast 
with the vast majority of modern work on the subject. 
However, the validity of our simplified model of cos- 
mic ray contribution requires further tests and compari- 



son with real supernova remnants and more sophisticated 
schemes modeling them. 

The value of the energy density found for relativis- 
tic particles Eqr can be used to estimate the intensity of 
magnetic field in the remnant. Assuming that Eqr and 
B^/^TT should be of the same order of magnitude, we de- 
rive B = 500+200 /iC This value appears to be somewhat 
high for this SNR, but still within reasonable limit s. For 
examp le, using the correlation pcsmVj k B^ from IVinkI 
( |2006[ Fig. 6) one might expect B < 100 juG. 

The list of the typical properties of these models is 
also presented in Table |2] In this case the difference 
between electron and ion temperatures is not so dras- 
tic: Te/Tj ^ . 03, w hich is still in agreement with 
^Gh avamian et al.l ( l2000l) . where at y = 4/3 (relativistic 
equation of state) the SNR plasma velocities (of ~ 4500 
km/s) correspond to T^ /T, ^ 10"^. 



5. Discussions and concluding remarks 

The analysis of the XMM-Newton X-ray spectrum of the 
SNR 0509-67.5 allowed us to estimate basic parameters 
of the object and reveals several peculiar features. 

First of all we obtained reasonable and consistent es- 
timates for the CSM density, using various methods: the 
EPIC MOS data fitting gives the estimates of the cir- 
cumstellar medium density of ncsM ^ 0.6 cm~^, the 
RGS spectral fitting of nitrogen abundance (which is 
not a product of thermonuclear explosion) yield «csm = 
0.4 - 0.8 cm"^. In the hydrodynamical simulations we 



Kosenko et al.: The X-ray Spectra Of The SNR 0509-67.5 9 

Table 2. Observed properties of the SNR 0509-67.5 and the properties of the numerical models: the one with the 
"basic" set up and one with a relativistic component (EqlTJ. The table columns are: R - visible remnant's radius (in the 
models this is the position of the contact discontinuity), t - the age, ncsM - circumstellar number density, v - velocity 
of the shocked ejecta (maximum value for the numerical models), kT^ - electron temperature (for numerical models it 
is the maximum temperature), kTj - maximum ion temperature, Mpe - an estirnate of the swept up iron mass (in solar 
masses), Nh - column density on the line of sight (derived from xspec. lArnaudll 1 9961 spectral fitting). 



Parameters 


R, pc 


t, yrs 


"CSM 


cm ^ 


V, km/s 


fcr„ keV 


kT„ keV 


MFe,Mo 


Nh 10''cm-2 


SNR 0509-67.5 


3.6 


<500 


0.4- 


0.6 


4900 ± 400 


2.5-3.6 


- 


-0.05 


- 


W7(basic) 


3.6 


400 


0.1 




<4600 


0.5 - 3.6 


150 - 2500 


-0.12 


0.7 + 0.3 


c050403m(basic) 


3.7 


400 


0.1 




<4200 


1- 1.5 


30-40 


-0.36 


5.6 + 0.3 


W7(CR) 


3.6 


400 


0.1 




<4300 


1.8-1.9 


20-36 


-0.12 


2.1+0.3 


c050403m(CR) 


3.8 


400 


0.1 




<4700 


2-45 


30 - 300 


- 0.36 


2.6 ± 0.3 



have found similar, but somewhat lower value of «csm = 
0.1 cm^l 

The abundances deduced from a single ionization 
timescale NEI (Fig[3]l show that swept up amount of iron 
in the remnant is low. Moreover the model was unable 
to reproduce the observed Fe K line, so the line was fit- 
ted separately. The centroid of Fe K line is located at 
6.50 + 0.05 keV, so it should be produced by low ionized 
(< XVII) ions of iron. This could mean that the reverse 
shock of the remnant just recently reached and heated up 
the iron core of the supernova. The Fe K feature enables 
us to make a crude estimate of the amount of heated iron 

ofMpe -O.IMq. 

In a framework of a NEI model with fixed single 
ionization timescale, abundances of the species lighter 
than aluminum tend to be overestimated, meanwhile, the 
amount of he avy species ( i ron g roup) tends to be un- 
derestimated dHughes et al.l 119981) . For example, in the 
W7 model the Mg/Si line flux ratio is higher than ob- 
served (Fig|6] Fig|8]l, while Figl3] shows that the Mg/Si 
abundance ratio is lower compared to the observed one. 
The sam e is true f or the iron abundance: iBadenes et alj 
(2008); Rest et al.l ( l2008h showed that the remnant has a 
very large amount of iron, while the abundance of Fe, de- 
duced from this simple fitting is considerably low. Note, 
however, that probably most of the Fe has not yet been 
heated by the reverse shock. 

In this context (using only abundances derived from 
single ionization timescale NEI model) it is hard to de- 
termine the precise explosion mechanism. If we consider 
only the species from Si to Ar, then a delayed-detonation 
model with Mp e - 0.9 Mq and an explosion energy of 
1.4 X 10^' ergs JWooslev et al.ll2007l) seems more prefer- 
able (FigO. The same tendency is confirmed by compar- 
ison of the observed spectra and the spectra from numer- 
ical simulations, performed in this study. T his conclusion 
is in agreement wi t h the r esults obtained bv lBadenes et alj 
(l2008l) : lRestetal.l(l2008h . 



A separate study of the remnant's bright southwest re- 
gion shows that this spot appears to be more metal-rich: 
our fitting shows that it is at least twice more abundant in 
metals than the average metallicity throughout the rem- 
nant. This points to either an asymmetric explosion of 
the SNR or asymmetri c CSM. Meanwhile note, that in 
iGhavamian et al.l (l2007l) it was mentioned that this region 
has a blueshifted velocity excess. Which is a hint for an 
asymmetric explosion. 

The numerical models are in agreement with the over- 
all properties of SNR 0509-67.5 inferred from the studies 
presented in Section|3j which are as following: the CSM 
density pcsM - 3 x 10"-^ g/cm\ the age f = 350 - 400 
years, and the plasma velocities v ^ 4000 - 5000 km/s. 
The ratio of electron to ion temperatures of TdTj ^ 0.01 
is in satisfact ory agreement with the expansion velocity of 
the remnant dGhavamian et alj2000l) . but h eavily depends 
on th e assumed explosion model (see also iBadenes et al.l 
|2006|) . Ejecta abundances govern electron population n,. 
and thus, the X-ray spectral shape. The latter is used to 
find the "best-fit" value of TelTi. 

The low value of r^/r, together with almost unde- 
tectable power-law emission point to the yet ineff'ective 
role of plasma instabilities and/or magnetic field at this 
stage in this object. Nevertheless, the simulations also in- 
dicate that the energy of relativistic particles is probably 
not negligible and may contribute up to ~ 40% of thermal 
energy, so the magnetic field (provided that Eqr ~ B^/d,n) 
might be slightly higher, ~ 500 /iG, than the exp ected 
value of -■ 100/iG based on the relation reported bv lVinkI 
(|2006|) . but within reasonable limits. 

We have found satisfactory hydrodynamical models 
for the remnant, but many problems still exist. 

In our synthetic spectra some of the line centroids and 
line fluxes ratios do not perfectly fit the observations. This 
is partly due to limited set of explosion models available 
and to still appreciable number of free parameters (elec- 
tron temperature, thermal conduction, magnetic fields, 
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cosmic ray particles, and other non-thermal and relativis- 
tic effects). For instance, more energetic explosion model 
(£ > 1.4 X 10^') could produce a supernova remnant at 
the age of 350 - 400 years with required high ejecta ve- 
locity, observed radius, and more acceptable higher value 
of CSM density (up to lO"^'* g/cm"'). Nevertheless such a 
model should contain a lot of ^''Ni and small amount of Si 
and S. The X-ray spectra produced by such an explosion 
model disagree with the observed one which contains a 
very well pronounced Si K and S K lines. 

Moreover, one of the shortcomings of the code is that 
for the calculation of the X-ray emission, and for the re- 
production of the observed spectrum around 1 keV, we 
need to update the atomic data for Fe L code. A task that 
we plan to undertake in the near future. 

To create more or less reliable models for supernova 
remnants, further developments of the physics in hydro- 
dynamical simulations are required. Self-consistent treat- 
ment of cosmic -rays acceleration is necessary. 3D simula- 
tions may shed some light on the influence of various in- 
stabilities on the SNR dynamics. Further analysis of high 
resolution X-ray spectra is needed to evaluate the condi- 
tions in other SNRs as well. 
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